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We investigate charmonium production in Pb+Pb collisions at LHC beam energy £'iab=2.76 A TeV at fixed- 
target experiment ( ^fs^m=^2 GeV). In the frame of a transport approach including cold and hot nuclear matter 
effects on charmonium evolution, we focus on the anti-shadowing effect on the nuclear modification factors 
Raa and for the Jlj/ yield and transverse momentum. The yield is more suppressed at less forward rapidity 
(yiab -2) than that at very forward rapidity (viab -4) due to the shadowing and anti-shadowing in different 
rapidity bins. 
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I. INTRODUCTION 

Recently a fixed-target experiment using the LHC beams 
has been proposed iQl], where the study on quarkonia in nu¬ 
clear collisions becomes specifically important, due to the 
wide parton distributions in phase space which is helpful to re¬ 
veal the charmonium production mechanism 101. Correspond¬ 
ing to the LHC beam energy £’iab=2.76 A TeV, where A is the 
nucleon number of the incident nucleus, the center-of-mass 
energy sfs^-12 GeV is in between the SPS and RHIC en¬ 
ergies, and a quark-gluon plasma is expected to be created in 
the early stage of heavy ion collisions. Taking into account 
the advantage of high luminosity in fixed-target experiments, 
which is helpful for detailed study of rare particles, the 
yield in Pb+Pb collisions at £’iab=2.76 A TeV per LHC run 
year is about 100 times larger than the T/i/^ yield in Au+Au 
collisions at ViArv=62.4 GeV per RHIC run year H. With 
the high statistics, one may precisely distinguish between dif¬ 
ferent cold and hot nuclear matter effects on charmoniunytro- 
duction |I3]. As is well known, the shadowing effect Hilt], 
namely the difference between the parton distributions in a 
nucleus and in a free nucleon, depends strongly on the par- 
ton momentum fraction x. Since x runs in a wide region, 
0.001 < X < 0.5, in the fixed-target experiments, it provides 
a chance to see clearly the shadowing effect on the charmo¬ 
nium distributions in different rapidity bins. In this paper, we 
study the shadowing effect on the nuclear modification factors 
for yield and transverse momentum in Pb+Pb collisions 
at LHC beam energy £’iab=2.76 A TeV. 

II. EVOLUTION OF QUARK-GLUON PLASMA 

The medium created in heavy ion collision at sJsnn - 72 
GeV is assumed to reach local equilibrium at a proper time 
To =0.6 fm/c i^, its consequent space-time evolution can be 
controlled by the ideal hydrodynamic equations, 

= 0 , 

5 ;./ = 0 , ( 1 ) 

where Tf,y = (e + pjUf^Uv - gfivP, j/i - nu^, m^, e, p and n 
are respectively the energy-momentum tensor, baryon current. 


four-velocity of the fluid cell, energy density, pressure and 
baryon density of the system. The solution of the hydrody¬ 
namic equations provides the local temperature T (x), baryon 
chemical potential p{x) and fluid velocity u^ix) of the medium 
which will be used in the calculation of the charmonium sup¬ 
pression and regeneration rates |01. Taking the assumption of 
Hubble-like expansion and initial boost invariance along the 
colliding direction for high energy nuclear collisions, we can 
employ the well tested 2+1 dimensional version of the hydro¬ 
dynamics in describing the evolution of the medium created at 
xjsNN = 72 GeV. Introducing the proper time t = and 

space-time rapidity rj = 1/2 In [(f 3- z)l{t - z)] instead of the 
time t and longitudinal coordinate z, the conservation equa¬ 
tions can be simplified as |0] 

drE + VM = -{E + p)/t, 

drM^ + V(M;cV) = -M^It - dxP, 

djMy + V(MyV) = -My It - dyP, 

drR + V(Rv) = -R/t (2) 

with the definitions E - (e 3- p)y^ - p, M = (e + p)y^\ and 
R - yn, where v and y are the three-velocity of the fluid cell 
and Lorentz factor in the transverse plane. 

To close the hydrodynamical equations one needs to know 
the equation of state of the medium. From recent studies on 
particle elliptic flow and shear viscosity, the matter created in 
heavy ion collisions at RHIC and LHC energies is very close 
to a perfect fluid Jgt]. Considering that the momentum inte¬ 
grated particle yield, especially for heavy quarkonia, is not 
sensitive to the equation of state, we follow Ref. M where 
the deconfined phase at high temperature is an ideal gas of 
gluons and massless u and d quarks plus 150 MeV massed s 
quarks, and the hadron phase at low temperature is an ideal 
gas of all known hadrons and resonances with mass up to 2 
GeVlfnl]. There is a first order phase transition between these 
two phases. In the mixed phase, the Maxwell construction is 
used. The mean field repulsion parameter and the bag con¬ 
stant are chosen as /if=450 MeV fm^ and B^^^-236 MeV to 
obtain the critical temperature Tc - 165 MeV ifl^ at vanish¬ 
ing baryon number density. Note that, when one calculates 
the rapidity or transverse momentum distribution of quarko¬ 
nia, the choice of the equation of state may result in sizeable 
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difference. 

The initialization of the hot medium is taken as the same 
treatment in Ref. Jst]. We use the final charged multiplicity 
to determine the initial entropy density. For yiNN=72 GeV, 
the charged multiplicity at central rapidity in center-of-mass 
frame is estimated to be dNchldrj = 515 based on the empiri¬ 
cal formula d: 

dNch , _ 

= 312.5 logio V^-64.8. (3) 

drj 

The initial baryon density is obtained by adjusting the en¬ 
tropy per baryon to be 250 Q. From the empirical relation 
o-NN = 29.797 + 0.141(ln iSl between the inelas¬ 

tic nucleon-nucleon cross section cr;viv in unit of mb and the 
colliding energy in unit of GeV, we have (j^n - 36 

mb at yfst^-12 GeV. These initial conditions lead to a max¬ 
imum medium temperature 7’o=310 MeV at the initial time 
to=0.6 fm/c. The medium maintains local chemical and ther¬ 
mal equilibrium during the evolution. If we do not consider 
the charmonium interaction with the hadron gas, the charmo- 
nium distributions in the final state will be fixed at time 
corresponding to the critical temperature Tc of the deconfine¬ 
ment phase transition. 


III. CHARMONIUM TRANSPORT IN QUARK-GLUON 
PLASMA 


Since a charmonium is so heavy, its equilibrium with the 
medium can hardly be reached, we use a Boltzmann trans¬ 
port equation to describe its phase space distribution function 
f'vix, p|b) in heavy ion collisions at impact parameter b. 


+ D^i, (4) 


where the loss and gain terms C^(x, p|b) and D^j(x, p|b) come 
from the charmonium dissociation and regeneration in the cre¬ 
ated hot medium. We have neglected here the elastic scat¬ 
tering, since the charmonium mass is much larger than the 
typical medium temperature. Considering that the feed-down 
from the excited states ij/' and Xc to the ground state Jlip S 
happens after the medium evolution, we should take transport 
equations for 'P = 7/i/r, if/' andxc when we calculate the Jfif/ 
distribution /y/,^ in the final state. 

Introducing the momentum rapidity y - 
1/2 In [(£-H P 2 )/(£’- Pj)] and transverse energy 

Ef - - p2 to replace the longitudinal momentum 

Pj and energy E = ^irp- + p“, the transport equation can be 
rewritten as 


sinh(y-77) 

cosh(j - r])dj H--h v, • V, 


h 


-crvy/j- -H 

(5) 


with the dissociation and regeneration rates p|b) = 

C^i(x, plbj/Zi, and p^i(x, p|b) = D^i(x, plbj/F,, where the third 
term in the square bracket arises from the free streaming of 
T* with transverse velocity v, =_p, /Et which leads to a strong 
leakage effect at SPS energy lU^ . 


Considering the gluon dissociation T' + g c-Hcin the 
quark-gluon plasma, the dissociation rate a can be expressed 
as 

(6) 

where Eg is the gluon energy, Fg^i - ^Jipkp - = pk 

the flux factor, and fg the gluon thermal distribution as a func¬ 
tion of the local temperature ^(.rlb) and fluid velocity M^(x|b) 
determined by the hydrodynamics. The dissociation cross sec¬ 
tion in vacuum cr^ip/p, k, 0) can be derived through the oper¬ 
ator production expansion (OPE) method with a perturbative 
Coulomb wave function lUTI - l^ . However, the method is no 
longer valid for loosely bound states at high temperature. To 
reasonably describe the temperature dependence of the cross 
section, we take the geometric relation between the averaged 
charmonium size and the cross section, 

(r^hiT) 

(Tg^ip, k, T) = l^^cr,4-(p, k, 0). (7) 

The averaged radial square {r~)^{T) is calculated via potential 
model iItIT with lattice simulated heavy quark potential ii 
at finite temperature. When T approaches to the charmonium 
dissociation temperature Ty, the averaged radius square and in 
turn the cross section go to infinity, which means a complete 
charmonium melting induced by color screening 0. Using 
the internal energy U as the heavy quark potential V, the dis¬ 
sociation temperautre Ty is calculated to be 2.17’c, \ .\6Tc and 
1.127’c for Jlij/,Xc and i(/', respectively 10 . 

The regeneration rate p is connected to the dissociation 
rate a via the detailed balance between the gluon dissoci¬ 
ation process and its inverse process |2^ . To obtain 
the regeneration rate, we also need the charm quark distri¬ 
bution function in medium. Although the initially produced 
charm quarks would carry high transverse momentum, they 
lose energy (momentum) when passing through the medium. 
Considering the experimentally observed large op en charm 
quench factor 11261 - 1^ and elliptic flow 11^ 1^ . we take 
as a first approximation a kinetically thermalized momentum 
spectrum for the charm quark distribution /c(x, q|b). Neglect¬ 
ing the creation and annihilation of charm-anticharm pairs in¬ 
side the medium, the spacial density of charm quark number 
Pc(-r|b) = f d^qlil/dy’fcix, q|b) satisfies the conservation law 

d,ipcu^)^0 ( 8 ) 

with the initial density determined by the nuclear geome¬ 
try pc(xo|b) = TAix,)TB(x, - b) cosh 7]/Toda^^/dr], where 
PA.sif^dz are the thickness functions, and 
dcr^^ I dr] is the charm quark rapidity distribution in p+p col¬ 
lisions. 

For the regeneration rate p, we also considered the canon¬ 
ical effect which is shown to be important in explaining the 
suppression of strange mesons 0. When there are only few 
pairs or even less than one pair of charm quarks produced in 
an event, one need to consider the canonical effect to guaran¬ 
tee the exact charm number conservation. Taking into account 
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the fact that the charm and anti-charm quarks inside a pair are 
produced at the same rapidity, we simply multiply the regener¬ 
ation rate g in a unit rapidity bin by a canonical enhancement 
factor ll^ 


Cc-c^l + llidNaldy). (9) 

To take into account the relativistic effect on the dissocia¬ 
tion cross section to avoid the divergence in the regeneration 
cross section, we should replace the charmonium binding en- 
er gy b y the gluon threshold energy in the calculations of a and 

gol. 

In the hadron phase of the fireball with temperature T < Tc, 
there are many effective models that can be used to cal¬ 
culate the inelastic cross sections between charmonia and 
hadrons 0. For 7/i/r the dissociation cross section is about 
a few mb which is comparable with the gluon dissociation 
cross section. However, considering that the hadron phase ap¬ 
pears in the later evolution of the fireball, the ingredient den¬ 
sity of the system is much more dilute in comparison with the 
early hot and dense period ||3l. Taking, for instance, the re¬ 
generation processes c + c g + Jftfr in quark matter and 
D + D* —» n + Jlxjj in hadron matter, the density ratio be¬ 
tween charm quarks at initial temperature Tq - 310 MeV and 
D mesons at critical temperature Tc - 165 MeV is around 
30. Considering further the life time of the quark matter ~ 6 
fm/c and the life time of the hadron matter ~ 2 fm/c calculated 
from the hydrodynamics in Section|II] we neglect the charmo¬ 
nium production and suppression in hadron gas, to simplify 
the numerical calculations. Note that, the suppression and re¬ 
generation in hadron gas may become important for excited 
charmonium states 11^ . 

The transport Muation can be solved analytically with the 
explicit solution Jj, 1^ 

h (Pr, y, X,, T], t) = fy (p,, y, X,(to), //(to). To) 

Jto A(t') 

X X' Q'T(P„.v,x,(T"),//(T"),r") ^ 


in the center of mass frame is /( y^NN /2) ~ 0.35 fm/c, 

which is compatible with the charmonium formation time but 
shorter than the QGP formation time to = 0.6 fm. Therefore, 
all the cold nuclear matter effects can be reflected in the initial 
charmonium distribution fp at time to. We take into account 
nuclear absorption, nuclear shadowing and Cronin effect. The 
initial distribution in the solution (fTOl) can be obtained from a 
superposition of p+p collisions, along with the modifications 
from these cold nuclear matter effects. 

The nuclear absorption is important in explaining the 7/i/' 
suppression in p+A and A+A collisions at low energies. It 
is due to the inelastic collision between the initially produced 
charmonia and the surrounding nucleons, and its effect on the 
charmonium surviving probability can be described by an ef¬ 
fective absorption cross section (Tabs- The value of (Tabs is usu¬ 
ally measured in p-i-A collisions and is several mb at SPS en¬ 
ergy. Since the nuclear absorption becomes weaker at higher 
colliding energy due to the shorter collision time llj^l^ . we 
take (Tabs=2 mb at £’iab=2.76 A TeV and the nuclear ab¬ 
sorption factor 

C -crabs(f”p(zA)*+/ifpfcXt-b)*) /n, 

5abs=e i-'' / ( 12 ) 

The Cronin effect broadens the momentum distribution of 
the initially produced charmonia in heavy ion collisions H. In 
p-i-A and A-t-A collisions, the incoming partons (both gluons 
and quarks) experience multiple scatterings with surrounding 
nucleons via soft gluon exchanges. The initial scatterings lead 
to an additional transverse momentum broade ning of partons 
which is then inherited by produced hadrons 11^ . Since the 
Cronin effect is caused by soft interactions, rigorous calcula¬ 
tions for the effect are not available. However, the effect is 
often treated as a random motion. Inspired from a random- 
walk picture, we take a Gaussian smearing SEIl] for the 
modified transverse momentum distribution 

_jsjM 1 n ~Pf~ 

/y (x,p,ZA,ZB|b) =- - I (f p;e“*"'/y (|Pf-p;i,Pz)5abs, 

Tiagi^l J 

(13) 

where 


with 

X,(t') = X, - Vt [TCOsh(y - t]) - t'A(t')] , 

//('!‘0 - y - arcsin (t/t' sinh(y - rj )), 

A(t') ^ yjI + (t/t')^ sinh^(y - 77). (11) 

The first and second terms on the right-hand side of the so¬ 
lution ([Tol l indicate the contributions from the initial produc¬ 
tion and continuous regeneration, respectively, and both suf¬ 
fer from the gluon dissociation in the medium. Since the re¬ 
generation happens in the deconfined phase, the regenerated 
quarkonia would have probability to be dissociated again by 
the surrounding gluons. The coordinate shifts x, —> X, and 
77 —> // in the solution (fTOl i reflect the leakage effect in the 
transverse and longitudinal directions. 

For fixed-target nuclear collisions at E]^\,-2.76 A TeV, the 
collision time for the two Pb nuclei to pass through each other 


l(x, za, Zij|b) = - ( r p(z, Xt)dz + f p(z, Xt - b)(/zj 
P \J~co Jzb / 

(14) 

is the path length of the two initial gluons in nuclei before 
fusing into a charmonium at x, za and zb, is the aver¬ 
aged charmonium transverse momentum square gained from 
the gluon scattering with a unit of length of nucleons, and 
/y^(P) is the charmonium momentum distribution in a free 
p-Hp collision. The Cronin parameter OgN is usually extracted 
from corresponding p-fA collisions. Considering the absence 
of p-i-A collision data at -s/Jm- 72 GeV, we take (2pM=0.085 
(GeV/c)^/fm from some empirical estimations As 

a comparison, for collisions at SPS (-y/sNN ~ 20 GeV) and 
RHIC (y^NN = 200 GeV) we take Og^ - 0.075 lit] and 
0.1 1 I 43 I] (GeV/c)^/fm, respectively. 

Assuming that the emitted gluon in the gluon fusion process 
g + g T' + gis soft in comparison with the initial gluons 
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and the produced charmonium and can be neglected in kine¬ 
matics, the charmonium production becomes a 2 —> 1 process 
approximately, and the longitudinal momentum fractions of 
the two initial gluons are calculated from the momentum con¬ 
servation. 


Xl,2 



e 


±y 


(15) 


The free distribution can be obtained by integrating 

the elementary partonic process, 

d(T^^ r dcr on —>¥0 

= J dygXiX2fg(XulUF)fg{X2,fiF) -, (16) 

where fg(x,iJ.F) is the gluon distribution in a free proton, yg 
the emitted gluon rapidity, dcrgg^^ig/dt the charmonium mo¬ 
mentum distribution produced from a gluon fusion process, 
and jjF the factorization scale of the fusion process. 

Now we consider the shadowing effect. The distribution 
function /,(x,yUf) for parton i in a nucleus differs from a su¬ 
perposition of the distribution f(x,fiF) in a free nucleon. The 
nuclear shadowing can be described by the modification fac¬ 
tor Rj - fil(Afj). To account for the spatial dependence of 
the shadowing in a finite nucleus, one assumes that the in¬ 
homogeneous shadowing is proportional to the parton path 
length through the nucleus 11 ^ . which amounts to consider 
the coherent interaction of the incident parton with all the tar¬ 
get partons along its path length. Therefore, we replace the 
homogeneous modification factor /?,(x,/i/r) by an inhomoge¬ 
neous one S 


Ri(x,fiF,^,) = 1 +A(Ri(x,fiF)-l)TA(xr)ITAB(0) (17) 

with the definition 7’.4B(b) = f d^x,TA{xt)T— b). We em¬ 
ploy in the following the EKS98 package ol to evaluate the 
homogeneous ratio R,, and the factorization scale is taken as 

Rf = yjml+pj. 

Replacing the free distribution fg in (fTbl) by the modified 
distribution = AfgRg and then taking into account the 
Cronin effect jl3b . we finally get the initial charmonium dis¬ 
tribution for the solution (fTOt . 

/y(xo, p|b) = r dZAdZBPA(tit,ZA)PB(t<‘r,ZB) 

E,Tq J 

X 'Rg(Xi,pF,X,)'Rg(X2,PF,^l - b) 

—NN 

X (x,p,ZA,ZB\b)Sabs- (18) 

Now the only thing left is the distribution in a free p-rp 
collision which can be fixed by experimental data or some 
model simulations. 


IV. NUMERICAL RESULTS 

The beam energy Eiab^ 2.76 A TeV in fixed target experi¬ 
ments corresponds to a colliding energy Vsnn= 72 GeV, and 


the rapidity in the center-of-mass frame is boosted in the lab¬ 
oratory frame with a rapidity shift Ay = tanh^'/lcms = 4.3. 
Let us first focus on the central rapidity region around ycms = 
0 in the center-of mass frame, which corresponds to yiab = 
4.3 in the laboratory frame. The centrality and momen¬ 
tum dependent anti-shadowing for initially produced char- 
monia is reflected in the inhomogeneous modification fac¬ 
tor Rg for gluons. The longitudinal momentum fractions 

are xi ,2 = + pj/ V'^nn ~ 0.05 for the two gluons, 

which is located at the strong anti-shadowing region by 
some parametrization of parton distribution shadowing like 
EKS98 H, EPS08 61] and EPS09 61. The anti-shadowing 
changes not only the gluon distribution but also the charm 
quark production cross section used in the regeneration. Eor 
the process g + g —> c + c, the anti-shadowing for gluons 
leads to an anti-shadowing factor ~ (Rg)^ for the cross sec¬ 
tion. Considering that in peripheral collisions the regenera¬ 
tion is weak and its contribution is not remarkably affected 
by the anti-shadowing, we take a centrality averaged anti¬ 
shadowing factor for the cross section to simplify the numer¬ 
ical calculation for regeneration. Estimated from the EKS98 
evolution il , we take a 20% enhancement of the charm quark 
production cross section compared to free pH-p collisions. 
Erom EONLL calculation 61 . the upper limit for dcr^^ jdy 
is 0.047 mb at yiNN=62.4 GeV. Note that the experimental 
data for charm quark cross section in free pH-p collisions are 
close to the upper limit of perturbative calculation, we take 
dcr^^ jdy - 0.05 mb at ysNN=72 GeV. After taking into ac¬ 
count the anti-shadowing effect in Ah-A collisions, it becomes 
0.06 mb. Eor p+p collisions, we assume a constant hidden to 
open charm ratio {dcr^i jdy) j {dcTcc jdy)-const at any colliding 
energy. Erom the ratio extracted from the RHIC data ifs^ . 
we have dcrjj^Idy=0.35 pb at yiNN=72 GeV. The transverse 
momentum distribution for Jjij/ in free p-t-p collisions can be 
simulated by PYTHIA [HI and the mean transverse momen¬ 
tum square is {p?)pp = 2.7 (GeV/c)^. 

Pig[T] shows our calculated centrality dependence of y/i/r 
nuclear modification factor Raa - j in Pb-i-Pb 

collisions at LHC beam energy £’iab=2.76 A TeV in labora¬ 
tory frame (GeV in center-of-mass frame) at for¬ 
ward rapidity yi^b = 4.3 (central rapidity ycms=0), where 
and are charmonium yields in p-i-p and A-i-A col¬ 
lisions, and Vcoii and Apart are numbers of binary collisions 
and participants. Eor comparison, we show also the RHIC 
data at ^]Jm — 62.4 GeV if^ at central rapidity. Since 
the shadowing/anti-shadowing effect is still an open question, 
and its degree depends strongly on the models we used, we 
show in Eigditwo calculations for the total Jlij/ Raa in Pb-t-Pb 
collisions at ysNN=72 GeV, one is with the above discussed 
anti-shadowing, and the other is without anti-shadowing. The 
hatched band is due to this uncertainty in the anti-shadowing. 
With increasing collision centrality, the initial contribution 
drops down, while the regeneration goes up. The canonical 
effect is important in peripheral collisions where the number 
of charm quark pairs is less than one and the inclusion of the 
canonical effect enhances sizeably the charmonium yield. In 
most central collisions, the regeneration can contribute about 
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FIG. 1. (Color online) The centrality dependence of the Jlijj nuclear 
modification factor at very forward rapidity yi^b = 4.3 (ycms=0) in 
Pb+Pb collisions at LHC beam energy E\a,=2.16 A TeV. The hatched 
band is the model result with the upper and lower borders corre¬ 
sponding to the calculations with and without anti-shadowing effect. 
The RHIC data are for Au-l-Au collisions at ycms=0. 


25% to the total charmonium yield. The anti-shadowing at 
very forward rapidity in the laboratory frame (central rapidity 
in the center-of-mass frame) enhances the charm quark cross 
section and in turn the initial charmonium yield by a factor of 
1.2. As a consequence, the enhancement factor for the regen¬ 
erated charmonium number is 1.2^ = 1.44 which leads to a 
strong charmonium enhancement! If we do not consider the 
anti-shadowing effect on the charmonium regeneration and 
initial production, the total Raa is significantly reduced. 


Figl2]we show the Jjif/ nuclear modification factor ifs^ 


rAA 


{P^AA 

^P^)pp 


(19) 


in Pb-t-Pb collisions at beam energy £’iab=2.76 A TeV, where 
{P^)aa and {pj)pp are averaged 7/i/r transverse momentum 
square in PbH-Pb and p-i-p collisions at very forward rapidity 
yiab=4.3. If we neglect the contribution from the regeneration 
and consider only the initial production, the ratio paa goes 
up monotonously with centrality due to the Cronin effect and 
leakage effect . The inclusion of regeneration (upper bor¬ 
der of the band) remarkably reduces the averaged transverse 
momentum, because the regenerated charmonia possess a soft 
momentum distribution induced by the charm quark energy 
loss. Since the degree of regeneration increases with central¬ 
ity, the increased soft component leads to a decreasing paa 
in most central collisions. The canonical effect can reduce the 
paa further, since it enhances the regeneration especially in pe¬ 
ripheral collisions. However, we should note that, the assump¬ 
tion of charm quark thermalization indicates a full energy loss 
and it may not be reached in peripheral and semi-central col¬ 
lisions at beam energy £’iab=2.76 A TeV. When we switch 
off the anti-shadowing (lower border of the band), both the 
hard component controlled by the initial production and the 
soft component dominated by the regeneration would be re¬ 
duced. Considering that the enhancement factor resulted from 
the anti-shadowing is 1.2 for the initial production but 1 .1? for 
the regeneration, the stronger anti-shadowing in the soft com¬ 
ponent leads to the slight difference between with and without 
considering the anti-shadowing, shown in Figl2] It is obvious 
that compared to the nuclear modification factor Raa for the 
yield, the modification factor paa for the transverse momen¬ 
tum is less sensitive to the shadowing effect ifs^ . 




FIG. 2. (Color online) The centrality dependence of the Jjifi nu¬ 
clear modification factor p^a at forward rapidity yiab=4.3 Cycm.s=0) in 
Pb-l-Pb collisions at LHC beam energy £’iab=2.76 A TeV. The upper 
and lower borders of the band correspond to the calculations with 
and without anti-shadowing effect. 

To see more clearly the charmonium production mecha¬ 
nism, we turn to the transverse momentum information. In 


FIG. 3. (Color online) The centrality dependence of the double ra¬ 
tios and for Jlijj yield and transverse 

momentum in Pb-l-Pb collisions at LHC beam energy £iab=2.76 A 
TeV. The upper and lower borders of the two bands correspond to the 
calculations with and without shadowing and anti-shadowing effects. 

From the simulations of parton distributions in cold nu¬ 
clear matter ifM WA. l4^ . the nuclear shadowing region is lo- 
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cated at very small x. In the following we consider the 
shadowing and see its difference from the anti-shadowing in 
Jjijj Raa and rAA in fixed-target Pb-i-Pb collisions. The max¬ 
imum 7/i/r rapidity in the center-of-mass frame is = 

cosh^^ [ V'^nn/ at y[s^-12 GeV. Considering 

the expected amount of measured events, we focus on the 
backward rapidity region around ycms = -2 which corre¬ 
sponds to the less forward rapidity yiab = Ay-i-ycms = 4.3-2 = 
2.3 in laboratory frame. From the kinematics, the momen¬ 
tum fractions for the two gluons involved in the gluon fu¬ 
sion process are xi = (+ pjI V'^nn)^^ = 0.35 and 

X 2 - ( Pf/ V'^nn)^'^ = 0.006. One is located in the 
EMC region and the other in the shadowing region 
leading to a reduction of 15% for the charm quark production 
cross section from EKS98 evolution ijst] (20% from EPS09 
NLO evolution JUI). Taking the same ratio of charm quark 
cross section between jcms = -2 and jcms = 0 calculated from 
EONLL ll^ and including the 15% shadowing reduction, we 
obtain c/cr^g^/c/y=0.01 mb at ycms=-2. Eor the medium evo¬ 
lution at this backward rapidity region, we initialize the en¬ 
tropy density to be half of that at central rapidity SHI which 
leads to a maximum temperature of 7'o=245 MeV. Eig|3]shows 
the two double ratios and of 

7/i/f, the upper and lower borders of the two bands correspond 
to the calculations with and without considering the nuclear 
shadowing and anti-shadowing. While the double ratio for 
the transverse momentum is not sensitive to the shadowing 
and anti-shadowing, as we discussed above, the strong anti¬ 
shadowing at yiab = 4.3 and shadowing at yi^b = 2.3 leads to a 
strong enhancement of the double ratio for the yield. Without 
considering the shadowing and anti-shadowing, the stronger 
charmonium suppression in the hotter medium at yiab = 4.3 
(ro=310 MeV) compared with the weaker suppression in the 
relatively colder medium at yiab = 2.3 (7’o=245 MeV) makes 
the double ratio less than unit. However, the inclusion of the 


yield enhancement due to the anti-shadowing at yiab = 4.3 
and the yield suppression due to the shadowing at yiab = 2.3 
changes significantly the behavior of the double ratio, it be¬ 
comes larger than unit and can reach 1.3 in most central colli¬ 
sions. Note that the rapidity dependent shadowing effect was 
used to qualitatively interpret the stronger suppression at for¬ 
ward r apidity than that at midrapidity in Auh-Au collisions at 
RHIC Hi HI. 


V. SUMMARY 

We investigated with a transport approach the charmonium 
production in fixed-target Pb-i-Pb collisions at LHC beam en¬ 
ergy £’iab=2.76 A TeV We focused on the rapidity dependent 
shadowing effect on the nuclear modification factors for the 
charmonium yield and transverse momentum. While the av¬ 
eraged transverse momentum is not sensitive to the shadowing 
effect, the anti-shadowing leads to a strong yield enhancement 
at very forward rapidity yi^b - 4, and the shadowing results in 
a strong yield suppression at less forward rapidity yiab -2. The 
double ratio between the nuclear modification factors Raa in 
the two rapidity regions amplifies the shadowing effect, it is 
larger than unit and can reach 1.3 in most central collisions. 

Erom the model studies on gluon distribution in nuclei, see 
for instance Refs. 1 H. I 4 U 4 I] . there are large uncertainties in 
the domain of large x (> 0.1), which is probably due to the 
unknown EMC effect. Erom our calculation here, the double 
ratio of the nuclear modification factor for 7/i/r yield is very 
sensitive to the gluon shadowing effect in different x region. 
A precise measurement of the ratio may provide a sensitive 
probe to the gluon distribution. 
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